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ABSTRACT 


A method utilizing vector representation is investigated for deter- 
mining a faulty element in passive and active networks by simple 
external measurements. 

A large system may be considered as an interconnection of a number 
of subnetworks. By utilizing the relationships between the magnitudes 
of a transfer function at various frequencies and the deviations of a 
circuit element, the fault simulation curves can be drawn. The fault 
identification regions are defined from the fault simulation curves. 

A fault identification matrix is menetructed corresponding to the de- 
fined fault identification regions. 

The fault identification matrix, when premultiplied by a vector 
whose components are measured from a network, yields another vector 
whose components identify a network element which is faulty. 

A test procedure for the fault identification method is presented 


and verified by experiments. 
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I. INTRODUCTION 


Recent advances in the design of electronic ere have resulted in 
an ever-increasing complexity in such systems. Despite the advances in 
component reliability and preventive maintenance, it appeared that a 
great many unpredictable or random failures do occur. Such complexity 
in modern equipments used by the industry and the military has resulted 
in a serious problem of providing adequate maintenance. One of the most 
difficult tasks is the location of malfunctions which, until recently, 
has been left entirely in the hands of technicians. 

Some techniques for fault identification have been studied by a 
number of investigators i 2 Se Gish On and a computer program for 
fault simulation of passive linear networks without mutual inductance 
has been given-in Ref. [2]. 

In this paper, the problem of fault identification, which requires 
the least number of external test points of the system, is investigated 
and a procedure which will automatically identify a fault or faults in 
passive and active linear networks is presented. Toward achieving this 
ultimate goal, the sensitivity of a network function due to element devi- 
ations is studied. The method, which requires matrix operations, is 
introduced. In this method, the fault identification regions at a set 
of selected test frequencies are defined depending on the network sensi- 
tivity and the fault identification matrix is generated. A vector is 
then obtained by premultiplying the fault identification matrix by 
another vector, and the resulting vector identifies the fault. 

The single element fault identification method is discussed and the 
combined element fault identification method is introduced by utilizing 


the single element fault identification method. 





The computer programs, which identify a fault element in a properly 
partitioned network, using the proposed method are written in FORTRAN IV 
and included in Appendix A. The implementation of the computer programs 


1S annexed in Appendix B. 





TJ. PRELIMINARY CONSIDERATIONS 


A. DEFINITION OF A VECTOR 
Consider a two-port network with the external terminals and variables 


as defined in Fig. 2-1]. The forward voltage ratio transfer function, 





2-port network. 


T12, 1s defined as the ratio of the output voltage to the input voltaye 


with the output current set to zero. 


A vee (2-1) 


The magnitude of 112, designated by |T12| , is given by either a 
function of the frequency with fixed element values or a function of an 
element value with a fixed frequency. That is, at<a set of n frequen- 
cies, with fixed element values, a set of n different magnitudes of a 
transfer function, T12, can be obtained. Similarly, corresponding to a 
set of k deviations of an element value, with a fixed frequency, a set 
of k different magnitudes of T1l2 can be determined. 

Let the network contain only the following types of passive elements: 
resistors, inductors and capacitors. Consider an element which has a 
nominal element value E with a deviation designated by aE. , where the 


subscript j implies the j-th deviation. 





With a change in element value AE, of E, n values of |112] can 
be obtained corresponding to n different frequencies. These values of 


|T12| are designated by 
[M12]. M12] «5 wore 4 M12] 5, eo, IT12] 


where the first subscript indicates the j-th deviation and the second 
Subscript indicates the i-th frequency. Next, with a nominal element 
value of E, n values of -|T12| can be obtained corresponding to n dif- 


ferent frequencies and designated by 
IT12] 9, : M2] 95 , 7" ITZ] 9 a ee IN2]| 5, 
The normalized value of M12]. is defined as 


712]. 
IT12] ). 
where i= 1, 2, --- , Nn. 
For example, with a deviation Ag. of E anda set of 5 frequencies, the 


5 values of the normalized |T12],. are 


ile alee IT12|,. 
ieee Yel pai eile 


With another deviation, AE, , of E and the same set of 5 frequencies, 


2 
the 5 values of the normalized IT12|,, are 


IT12|,. ITM12| 55 | I712| 56 
Maree zig, °° 2 600” > TMZ] 55 





Let a vector X be denoted by 


X = ( me X93 oar 5 Xa» ae 9 STE C ) (2-2) 


where X is an array of the n normalized values of {T12] with respect 
to a set ni ell (with a deviation of aE.) dhig ameons tan tac: 
ia ee. 
x. denotes ae (i = 1,2, --- , n) andc is unity; 
1 [T12] 9. 
X is n dimensional because the last component is always unity; 
the superscript (') means a transpose of an array of X or a matrix. 
B. CONSTRUCTION OF FAULT IDENTIFICATION MATRIX AND 
GENERATION OF VECTOR-Z 
Consider a simple two-port network that contains a resistor, R, and 


a capacitor, C, as shown in Fig. 2-2. Let the nominal value of R be 


1 ohm and that of C be 1 farad. The value of C can be changed from 





RiGee 252. 
2-port RC network. 


zero to infinity and 3 nonzero finite values of C are indicated as shown 


in Fig. 2-3. 


O20 U0 ee leo 





In general, consider that the value of a parameter can vary from 
zero to infinity and this range of values can be divided into a number 
of subdivisions, each of which shall be called the j-th deviation of the 
Berameter Value, wnere j = 1, 2, --- . For example, the value of C 
(between 0 to ~) is divided into four subdivisons, that is, the 1]-st 
Meyiation is 0 < € < 0.5 , the 2-nd deviation is 0.5 <C< 1.2 and 
sarily, J.2c =< C< 1.5 and 1.5 < C represent, respectively, the 
3-rd and 4-th deviation. - 


The normalized |T12] jis given by 


Mtg vw (223) 
[112] 95 |i + (we)? 


The values of the normalized |T12] corresponding to the 1-st deviation 
through the 4-th deviation are shown in Fig. 2-4 according to a set of 
3 test frequencies. The values of the normalized |T12| corresponding 
to the j-th deviation where j = 1,2,3,4 are best seen by a coordinate 
plane whose ordinate represents a Rat cHey axis and abscissa represents 
a normalized |T1l2| axis as shown in Fig. 2-5. | 

Let a "region-j", designated by R-j, represent an area on the fre- 
quency-normalized |T12| coordinate plane which corresponds to the j-th 
deviation of an element. For example, R-1 represents the area to the 
right of line A-O (Fig. 2-5) on the frequency-normalized |T12| coordi- 
nate plane and corresponds to the |-st deviation of C. R-2 represents 
an area between line B-O and line A-0, including line A-0, on the fre- 
quency-normalized |112] coordinate plane and corresponds to the 2-nd 


deviation of C as shown in the figure. 
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T-st 2-nd 3-rd 4-th 
deviation deviation deviation deviation 


Oi. 0.20 We eS 
















- 0.906 Oey o 
=) [ 7 785 5 
0.861 0.707 
w=2 ee ete ree % 
2 u ¢ 
w=4 X 0.845 Xa 0.68 X. 
Fig. 2-4. 
The normalized |712| corresponding to the j-th deviation. 
Wi 
6 B A 
; R-3 
(3-rd dev. ) 
, 
4 cas oa *3 a 
3 . R-4 R-2 R-] 
(4-th dev.) (]-st dev.) 
x 
2 a SS c i Xp 
x ‘] 
] <——— *] <— - 
0 
0 1.0 2.0 


"normalized |T12| 


dic 3. 


The frequency-normalized |712| coordinate plane. 
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A “normal region’ is defined as a region on the frequency-normalized 
[T12| coordinate plane where a corresponding element deviation is ina 
predefined normal range of deviation. The region outside of this normal 
region is defined as a "faulty region". For example, if the 2-nd devi- 
ation is a predefined normal range of deviation, R-2 is a normal region. 
The regions outside of R-2, that is, R-1, R-3 and R-4 are considered to 
be faulty regions. 

Let X. be a variable in a region R-j of domain x. If the value of 
Xs is in a normal region, the corresponding element is considered to be 
normal. For example, let R-2 be a normal region. If the values of x. 
where i = 1,2,3 is in R-2, the element C is’ considered to be normal. 
Otherwise the element with the values represented by X is considered 
to be faulty. . | 

Let a "boundary point value", abbreviated by B.P.V., represent a 
value of the normalized |T12| which lies between two regions. The 
boundary point values of a region may or may not belong to that region. 
For example, the boundary point values between R-1 and R-2 lie on line 
A-O and these values do not belong to R-1} but belong to R-2 as shown 
in Fig. 2-5. The 3 boundary point values corresponding to Ky» Xo and 
Xo in R-1 are 1.265 , 1.582 and 1.845 and belong to R-2. 

A “fault identification matrix", designated by M, is defined as a 
Matrix whose number of columns represent the number of regions and the 
number of rows represent the number of components of X as shown in 
fle 2-6. 


A vector '" Z" is defined as 


ee a (2-4) 





(No. of regions = No. of columns) 


Pelee re Ray ee- Rey 
w= =fon enn dennneen-- ponann nan a) Va 
a eee oad 
aa , | 
M = _ 3 | i 
] j ! j 
j ] ] ( 
es oe 
nnn fennnndenn------ jo-------- geen) Deg 


(No. of X components = No. of rows) 


Elgin 2-0. 


Fault identification matrix. 


A method for the construction of M and the generation of Z in the 


n-dimensional domain x with r-regions is presented as follows. The 


r-regions R-1, R-2, ----- , R-j), ----- » R-r and p boundary point values 
by» b,, ----- : b,, » assigned by positive numbers, are defined in Fig. 
2-7. 


Corresponding to the boundary point values, the r-regions can be 
arranged with respect to the vector component x. as’ shown in Table 2-1. 
Let the regions, arranged as shown in Table 2-1, be called the "identi- 
fication regions". 

A fault identification matrix, M, will be constructed corresponding 
to the defined identification regions. For every B.P.V. associated 
with each vector component, xX 5 a sub-fault identification matrix is 
formed. These sub-fault identification matrices, designated by My > Ms 


---- , have the same dimension as M. 


iis 
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R-1 
R-2 


Seay 


n-dimensional domain x with r-regions 





Table 2-1]. 
Identification regions. 
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The following rules are given for the construction of a fault identi- 
fication matrix. 
Rule 1: 
To construct a sub-fault identification matrix for a particular 
Xs the following procedures are used. 

(a). If the value of x in R-j is greater than the selected 
B.P.V., that is associated with Xa put +] in the i-th row of j-th column 
of the sub-fault identification matrix. 

(b). If the value of x in R-j is less than the selected 
Ber.¥., that is associated with Xe put -] in the i-th row of j-th column 
of the sub-fault identification matrix. 

(c). To generate the entry of the j-th column of the last 
row, multiply the selected B.P.V. by the entry value of the j-th colum 
of the i-th row and change its sign. 

(d). All entries other than those in the i-th row and the 
last row are set to zero. 

Rule 2: 
The fault identification matrix, M, is constructed as a matrix 
sum of all-sub-fault identification matrices. 
Using rules 1 and 2, the sub-fault identification matrices and 
M for the identification regions in Table 2-1 are constructed as follows. 


With respect to Xy> the sub-fault identification matrices are: 





for B.P.V. of D(p-1)? 


and for B.P.V. of Do» 


With respect to Xo» the sub-fault identification matrices are: 


for B.P.V.- of bo» 


stat) 


R-2 


ae pE |) 


R-2 


+b 


+b 
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oa Re] 


n+] 


n+] 


n+] 


a 





ena tor B.P.V. of b 


Finally, with respect to Xa? the sub-fault identification matrix for 


ferey. of b_ is 
mer p 


where Me means the 


The resulting M is 


M 


ie 
R-1 R-2 --- Rj --- R-r 
= 0 0. 0 
+] +] +] 
0 OF 0 
“by -b, -b, 


R-| R-2 --- R-j --- R-r 
= 0 0. 0 
0 0. 0 
-] =e +] 
+b +b -b 
p 
last sub-fault identification matrix. 
given by 
ee eueice Uy + Me 
R-| R-2 ~-- R-j --- R-r 
= Oy =e ee 
0. 0. ce 
-| -1. +] 
Ay Ao Ay. 
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n+] 


n+] 


nt+1 


a 





- b, + --- + b 


where A, = bip-1) ee eae t Bb > A 9 - by ; 


> = ire) + 2b 


and ae = oat) - eb. - b. + --- - b ; 
Next, the application of M for the generation of Z is discussed. 
Let the normalized values of |1T12| of a system be eee == 5 oN 


Des in the form of 


><| 


[ Ve oO aie east CON mr Ds 5 ie ] 
Z' is generated as follows. 


ZT 


i 
>< 
= 
am 
~ 

— 
*< 
nN 

I 

‘ 

I 

I 

I 

1 

1 

I 
< 
— 

es 
= 


I 
rc 
i 
bas 
= 
+. 
> 
—! 
nh 
Pas 
—! 
4" 
i 
m< 
s 
a 
_- 
nm 
i 
! 
M7 
~ 
he 
no 
> 


i ae 


= [ a 9 | Zo Soe tee. {5 ae 5 == s 9 Z J 


Z has the same number of components as the number of regions and 2 
serves as an indicator of the region in which the system is operating as 
given by the Fol lowing Rule: 
Rule 3: 

If a system is operating in R-j, the component 2, of Z is the alge- 
braically largest component in Z. 

The construction of M and the generation of Z are illustrated in 
the following examples. 


Example 2.1. Construct M and generate Z in the following one 


dimensional domain x with three regions. 


R-2 R-] R-3 
1 | 
0 1.0 
mre 5 .0 Se) Xa 10.0 50° < X 





Corresponding to the given boundary point values, the identification 


regions are as follows. 


R-1 R-2 R-3 
Xx Dre oe oe x50 10.0 < x 


Using rules 1 and 2, the sub-fault identification matrices and M for the 


identification regions are constructed. 


mreesub-fault identification matrix for 8.P.V. of 5.0 is 
R-| R-2 —  R=3 
my = +1, -l. : +1, 
-5.0 +5.0 -5.0 
and the sub-fault identification matrix for B.P.V. of 10.0 is 
R-] R-2 R-3 
Mm, = -l, -l. tale | 
+10.0 ar UGH 8,  =-10.0 
The resulting M is given by 
R-] R-2 R-3 
M = m, = mM, = 0. -2. ce | 
50 le 0 -15.0 


Suppose that the measured value of x is 13.0. X is in the form of 


Meee iox . 1. |} 
aoaOeel | 
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Thus, Z' is obtained 


Z' = X'M imeeO el. || 0. -2. a 
ceOmel5 0 15.0 


Mie es 


il 
rt 
nn 


[ 2 ; Zn» Z2 J 


in which 2. is the algebraically largest component. According to rule 3, 
the system is operating in R-3. The measured value. of x, indeed, belongs 


to R-3. 


If the measured value of x is 6.0, then 


Acie cro 1. I s ae 2. 
5.0 15.0 -15.0 


32 35 | 


it if 
ro co 
OT 


le as ae J 


The fact that zy = 5 , the largest of all the entries in Z', indicates 


that the system is operating in R-1. 
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Construct M and generate Z in the following three 


Example 2.2. 


dimensional domain x with four regions. 


pa ae ee ee ee ee x ea ee Terenas oe 
la ieee aA. 
+ 
~~ + NI 
x Pa * 
~~ — 
_ x x 
~ 
BE 9i hh Par eng | Sunroom emcee eno fa vet cme oma necro eee TL ane ec eR Qe as | rere ee aoe 
“I A f 
NI NI 
x ~~ 
= V 
OQ) oe re nn rn nr rm rr rrr rrr rrp 
A 
NI _ 
~*~ ~~ 
wo Joe eee oe) (a | ae 


R-2 
R-3 
R-4 


a 


The 


The boundary point values are 1.5 , 2.0 , 2.5 , 4.0 and 5.0. 


identification regions are defined as follows. 





2] 





With respect to Xy> the sub-fault identification matrices are: 


mor B.P.V. of 4.0, 


R-1 R-2 R-3 R-4 
m = [-l J a 
0 0. 0 0 
0 0. 0 0. 
4 4 -4 as 
eme nor B.P.V. of 2.0, 
R-1 22 re R-4 
My = i -] ] ] 
0. 0 0. 0 
0 0. 0 0 
= 2 a E 


With respect to Xo» the sub-fault identification matrices are: 


won b.P.V. of 2.5, 


R-1 R-2 R-3 R-4 
Mz = oF 0 0 0 
-1 “1, 1 1 
0 0 0. 0 

| 2.5 2.5 -2.5 2.5 
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ema for BoP. V. of 1.5, 


R-] R-2 R-3 
My, = OF 0. oF 
he ile. ik 
0. OF 0. 
oreo Sales -1.5 
With respect to Xo the sub-fault identification matrix 
5.0 is 
R-] R-2 R-3 
Me = ie 0. On 
OF 0. 0. 
-l. -l. -l. 
oi oe oe 
The resulting M is given by 
R-] R-2 
M = Mm, i M +m, +m + mp = 0. -2. 
OF OF 
=e -1. 
Me 12: 
If the measured components of X are X, = 3.0 5 X= 
X' becomes 
“eee co 8k Cl.) 
eee wis. co ce. 6. dM = CUdE 5. OS. 
[ Z1 > Zo > 
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R-4 
0; ] 
lw 2 
0. 3 
=1\. SE 4 
hore ber. 
R-4 
oR 1 
OF 2 
de 3 
-5. ip 4 
R-3 
2. 
2D. 
-1. 
ey 
O and Xo 
pom 1: ] 
23s Zy ] 


of 





The component Z, is the albebraically largest component. According to 
rule 3, the system is operating in R-1]. Similarly, if the measured 


components of X are X, = 3.0 , X, = 2.0 and x, = 6.0, Z is given by 


mee X Me=- [ 3. 2. 6. 1. )] M 


feces |. . 45. | 


Iezq Zo. Z, > Z, J 


The component Z4 indicates that the system is operating in R-4. 
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TIZ. JHE PROPOSED METHOD 


fee PURPOSE 

The purpose of the proposed method rot investigate a feasible method 
to determine whether fault detection and isolation can be accomplished by 
measuring only the magnitude of the system transfer frunction. In this 
method, the least number of test frequencies and test points, such as the 


normal input and output terminals of the system, will be used. 


B. BASIS FOR THE PROPOSED METHOD 
Consider again the 2-port network shown in Fig. 2-1. The open circuit 
impedance functions are represented by the equations of the voltages as 


linear functions of the currents in matrix form as follows. 


= (3-1) 


The network of the type being considered can be externally character- 
ized by a transfer function in terms of open circuit impedances. It is 
easily found that the forward voltage ratio transfer function, T1l2, and 
the reverse voltage ratio transfer function, T2l, (defined as the ratio 


of the input voltage to the output voltage with the input current set to 


zero) are 
Zz 
Tee A (3-2) 
11 
Z 
im a (53) 
22 
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In general, a transfer function can be expressed as a rational func- 


tion of the complex frequency, s, and a network element, E, as follows. 


(Ce (3-4) 


where zero, Ze and pole, P. , and K are functions of E and m and n are 
the number of zeroes and poles respectively. 


Taking the logarithm of Eq. (3-4), the following equation is obtained 


m n 
ineiiseey = In K+ » In(s+z.)- &§ In(s + p.) (3-5) 
i=] jel ’ 


And taking the derivative of Eq. (3-5), the following equation is 


obtained. 








d(In T(s,E)) = a ) 
m aver n Gore 
ees eee — E C(3-6) 
ounce i=] StZ. dE j=l s+ Ps dE 


A “transfer function sensitivity” is represented by Eq. (3-6). The 


transfer function sensitivity is closely related to pole and zero sensi- 


aides a2s 
tivities where aa and aa are defined as pole sensitivity and 


zero sensitivity respectively. In other words, the variation of T(s,—) 
may take the form of a change in location of a pole and/or a zero of 


Hise). 


“16 


a ee 
a - 


\ 





From Eq. (3-4) and Eq. (3-6), it can be seen that a change in any 
one of the pole-zero configurations will cause a change in one or possi- 
bly both of the voltage ratio transfer functions. Since T12(s,E) and 
T21(s,E) depend on the value of E, it is not hard to consider how a 
change of an internal element is related to a change of an external 
voltage ratio transfer function. 

This is the basic philosophy of the proposed method for fault 


identification. 


Meee PROPOSED TEST PROCEDURE 

A fundanental fault testing procedure for the proposed method is 
given as follows. 

Step ( Partition the system. 
Step ( Generate a symbolic system function. 
Step ( Select a set of test frequencies. 
Step ( Simulate a fault. 


Step ( Construct the fault identification matrix. 


Step ( Generate X. 


Generate Z. 


1) 
2) 
3) 
4) 
Step (5). Determine the fault identification regions. 
6) | 
7) 
Step (8) 
9) 


Step ( Identify the fault. 

(a). Single element. 

(b). Combined element. 

Step (2), step (3) and part of step (4) have been studied and a 
computer program of these steps for the passive linear network without 
mutual inductance was presented by J. D. Courville in Ref. [2]. The 


detailed discussion of these steps will be omitted. The other steps are 


discussed in detail in the next chapter illustrative examples. 


oy 





IV. NETWORK FAULT IDENTIFICATION 


A. NETWORK CONFIGURATIONS AND TRANSFER FUNCTIONS 

As illustrated in the foregoing discussions, there is a certain 
relationship between the transfer function and a network element devia-~ 
10N. 

To illustrate the proposed method, two examples for passive network 


configurations are presented. 


Example 4.1: 





1/Eal 


= Ry = 2 ohms 
}/Ea2 = L = 2 henry 
Ea3 = C = 3 farads 
1/Ea4 = Ro = 4 ohms 


Fig. 4-1(a). 
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Example 4.2: 


V 
V = 2 
] C cS L, 
ro Sz! 
1/Ebl = L, = | henry 
1/Eb2 = R = 1 ohm 
Eb3 = C = 1 farad 
1/Eb4 = Lo = | henry 
Fig. 4-1(b). 


For example 4.1, the transfer function is 


R 
Ti2 é 


2 
R,LCs + (R)RoC+L)s + Ry + Ro 


EPG Sse 5 


s 
= (4-1a) 
(st+0.262)(s+1.90) 


where Ky ise oes 


ao 





For example 4.2, the transfer function is 


Les 
lel é 


So 2 
L,L.Cs + RL, Cs + (Lj+L,)s + R 





eee (4-1b) 


(540.57) (s¢40.43s+1. 7623) 


where K, Ss lie 


The network semilogarithmic plots for the magnitude characteristic 
of the transfer function are shown in Fig. 4-2a for Eq. 4-la and Fig. 


4-2b for Eq. 4-1b. 


B. NECESSARY DATA 

All necessary data for the proposed method are based on a set of 
chosen test frequencies. 

A “fault simulation curve" is defined as a response curve of the 
normalized |T12| variations due to the simulated deviations of an 
element with respect to a chosen test frequency. The fault simulated 
data for |T12| can be obtained from the computer program in Ref. [2]. 
The fault simulated data for |T12] and the normalized |T12| are 
tabulated in Tables 4-la and 4-1lb for Eq. 4-la and Eq. 4-1b respectively. 

The fault simulation curves are drawn in Graph 4-la-1 thru Graph 


4-la-3 for example 4.1 and Graph 4-1b-1 thru Graph 4-1b-3 for example 4.2. 
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C. DETERMINATION OF FAULT IDENTIFICATION REGIONS 
The proper determination of the identification regions is the most 
laborious part of the proposed method. 


For example 4.1, a set of 3 test frequencies is chosen as follows. 


W) = 2.00 radians/sec 
Wo = ay radians/sec 
W = OO radians/sec 


For example 4.2, a set of 3 test frequencies is chosen as follows 


ee 3.20 radians/sec 
Wo = 0.70 radians/sec 
Note O205 radians/sec 


For a test frequency, there are eight basic identification regions, 
four normal regions and four faulty regions, because four elements are 
contained in the examples. That is, each element has two regions, a 
normal region and a faulty region. The four normal regions are not the 
concern of the propsed method. The four faulty regions may be sub- 
divided into many faulty regions to identify a fault. These subdivisions 
of the faulty regions depend on the transfer function sensitivity. The 
subdivided faulty regions are called the fault identification regions. 

The normal regions for the networks are defined as in Table 4-2a for 
example 4.1 and Table 4-2b for example 4.2 according to the fault simu- 
lation curves. The regions outside of the normal regions are considered 


to be faulty regions. The fault identification regions are defined as 
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Table 4-la. 


The magnitudes and the normalized magnitudes of 


T12 due to each element deviations for Example 4.1. 


20, 


= 


0.02 ) 


1.02 and 0.6647 at w 


( The nominal magnitudes are 0.05984 at 


0.1465 at w 
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Table 4-1b. 


The magnitudes and the normalized magnitudes of 


Tl2 due to each element deviations for Example 4.2. 


Creme 


0.05 ) 


0.7 and 0.04986 at w= 


( The nominal magnitudes are 0.11453 at w 
0.5964 at w 
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in Table 4-3a for example 4.1 and Table 4-3b for example 4.2. In Tables 
4-3a and 4-3b, an arrow "t" means "more than".and an arrow "+" means 


"less than". 


D. CONSTRUCTION OF FAULT IDENTIFICATION NATRIX 

The sub-fault identification matrices and the fault identification 
matrix, M, for the fault identification regions defined in Table 4-3a, 
are constructed step by step according to the rules given previously. 
The boundary point values for x where 1 = 1,2,3 are tabulated in Table 


4-4, 
ies 


Bree 


Table 4-4, 


The sub-fault identification matrices with respect to X, are repre- 
sented by mM M5» My and Wr , the sub-fault identification matrices with 
respect to x, are represented by Mos Me, Mz and Mo and the sub-fault 
identification matrices with respect to X_ are represented by Mg> Migs 
mM and M15: All of the sub-fault identification matrices and M are 


given as follows. 
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E. GENERATION OF X 

The components of X can be obtained by the measurements of the 
normalized |T12| corresponding to the chosen test frequencies. 

Suppose that it is desired to determine which element is faulty in 
a network accessible only at the input and output terminals. The 
diagnosis of the network malfunction can be accomplished using the fault 
identification matrix. 

The value of X will fall in one of the defined fault identification 


regions. 


F. GENERATION OF Z AND FAULT IDENTIFICATION 
The generation of Z and the fault identification are best seen by 
the following three cases for example 4.1. 
| Case 1. Let the measured values of the components of X at the 


selected test frequencies Wis Wo and Wo be 


X) = .6 
Xo = whi 
X3 = 9 


Z is generated as follows. 
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the components of X at 


lie 


ee 


-1.56 
is 


1.84 
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OO 


Case 3. If the measured values of the components of X at the 


selected test frequencies are 


Xy = 2 
X5 = lei 
Xa = 98 
fee = X'M 
= [| 2, 1.1 98 1. J M 


[ -4.34 4.34 -3.7 -4.3 3.18 4.78 


Pimms aiisch  =3.74 2.94 4.54 | 


7 EB 29 10 ah Za 


From the defined fault identification regions, Table 4-3a, the 
following facts are known for each case. 

Case 1. The algebraically largest component in Z is the component 2): 
The network is operating in R-l which means Eal is deviated less than 
-20 percent from the nominal value. In example 4.1, R, re out of toler- 
ance. 

Case 2. The component ze, in Z is the algebraically largest component. 
The network is operating in R-6. The element Ea2 is deviated more than 
60 percent from the nominal value. In example 4.1, L is out of tolerance. 

Case 3. The component Ze in Z is the algebraically largest component 
again but the noticeable fact is that the value of the component Z¢ is 
different from that of the previous case. The relative value of the 


component of Z may suggest the degree of violation of the faulty element. 


De 





The method discussed in the foregoing illustrations is called the 
Single element fault identification methcd. A combined element fault 
identification, abbreviated by C. E. F. I., method is a method in which 
more than one element fault is identified. 

The algorithm of the single element fault identification method is 
Suitable for computer implementation. If the fault identification 
regions are properly defined, then the construction of M and the genera- 
tion of Z are matrix operations. The fault identification is the 
comparison of the values of the components in Z. 

A program for the single element fault identification, designated by 
program |, is written as shown in Appendix A and its implementation is 
annexed in Appendix B. The results of the computer applications for 
example 4.1] are. shown in Fig. 4-3, Fig. 4-4 and Fig. 4-5. The same 
procedures are used for example 4.2 and the results of the computer 


applications are shown in Fig. 4-6, Fig. 4-/ and Fig. 4-8. 
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Fault identification regions 
printed by computer for Example 4.1. 
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V. COMBINED ELEMENT FAULT IDENTIFICATION 


Multi-element fault identification, abbreviated by M.E.F.1., is a 
somewhat complicated task. However, if mul ti-element fault simulation, 
abbreviated by M.E.F.S., for a certain system can be made, M.E.F.1. is 
possible. 

Consider the number of times of M.E.F.S. necessary for a M.E.F.1.. 
For a 3-element system, using the abbreviation for the combination of m 
elements out of n elements, nCm, the number of times to be simulated is 


as follows. 


3C1 = 3 --- one element fault combinations 
3C2. = See two elements fault combinations 
303 = |] --- three elements fault combination 


Including every element deviation of 10 times, each element deviation 
from -90 percent to +90 percent of the nominal value with 20 percent of 


Step size, the number of fault simulations necessary for a M.E.F.I. is 


aS many as 
sc 10 = 30 for one element fault simulations 
3C2 10 10 = 300 for two elements fault simulations 
3C3 10 10 10 = 1000 for three elements fault simulations 


1330 for total number of fault simulations 


Thus, a total of 1330 multi-element fault simulations are required even 
in a simple 3-element system. 
The number of possible M.E.F.S. required for a M.E.F.1I. is shown in 


Table 5-1 and it becomes astronomical as the number of elements increases. 
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No. of element NG Oiroiees tr. Oo. 
2 | 110 


3 | 1330 

4 14640 

5 161050 

6 1771560 
Table 5-1}. 


This astronomical increase of M.E.F.S. discourages one to make any 
further investigation. However, another method can take its place as is 
explained as follows. A complex system can be partitioned properly such 
that the probability of failure involving more than one element at a 
time is usually. small. 

In this section, the combined element fault identification (C.E.F.I.) 
method is introduced. A complex network shown in Fig. 5-1 is partitioned 
properly into two subnetworks shown in Fig. 5-2. Two subnetworks NW-1 
and NW-2 are regarded as independent and the single element fault identi- 


fication method can be applied to each network. 


Flaca |. 


6) 


* 








EE 


fie ae 


Suppose that the fault identification matrices for NW-1 and NW-2 are 


constructed and designated by 


M1 for Ni-1 
M2 for NW-2 


Let M be defined as 


M = y ! . 
ee (5=1)) 
0. M2 


The measured values of the components of X can be obtained with 


respect to a set of selected test frequencies for NW-1 and NW-2 independ- 


ently and designated as follows. 


Xyq2 Xyo2 Xqge cote > Xay> for NW-1 


Xo19 Xo9s Xog9 777mm > Xons for NW-2 


where n and m are the number of selected test frequencies for NW-1 and 


NW-2 respectively. Then Z are generated as follows. 
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—_——_— — = =o 


ee Me KY, 2" ke 7 ‘ = XI'M1, X2'M2 


O , M2 
= D2) 29 237m ag 2nd) “777 Zp 
where 

(ee Gee 

ee 2 3 777 Xi, |: J 
X2" = T Xy %XyQ X37 Xon t+ J 
Zi' = XI'M = [ 2, %y 24 oie ] 
Le = Xo M2 [ “(rdt1) =. Zh ] 


ris the total number of fault identification regions for both networks 
and rd and r-rd are the number of fault identification regions for NW-1 
and NW-2 respectively. Z will indicate two faulty elements. One from 
Z1 indicates a fault in NW-1 and the other one from Z2 indicates a fault 
in NW-2. 

This is the C.E.F.I. method for two elements fault identification. 
The C.E.F.I. method can be extended to isolate wil aelenant faults. 

For example, suppose that the network shown in Fig. 4-la is replaced 
by NW-1 and the network shown in Fig. 4-1b is replaced by NW-2 for 
C.E.F.I1. and the fault identification matrix for NW-1 is M1 and the 
fault identification matrix for NW-2 is M2. The resulting M is constructed 
by Eq. (5-1) and shown in Fig. 5-3. 


let the measured values of the components of X be 
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Xl' = [ X17 X02 X13 eee = wtOsoe 0.7 O.9 1. 

aa [ Xo) X09 X52 ete OG 294 50-8 Lael 

eee -emiOn6 07 0.9 |, 0.6 O.94 0.8 1. 
Z' becomes 


a= eM = XI'M1, X2'M2 


“we i 3.8 -3.18 2.94 2.7/4 -0.26 -1.46 3.14 -2.86 
2.98 MOC 0.96 =2.18 BiOd S207 3909 1.11 
SON S37 3.57. -0.83 -0.15 -0.39 1.45 -0.45 


2.67 | 
eulecyecan <3 -" ID Se Zon, J 
One fault identification region from R-] ----- R-12, represented by 
a Z45 respectively, indicates a fault in NW-1 and the other one 
fault identification region from R-13 ----- R-25, represented by 243 -- 


=" Zor respectively, indicates a fault in NW-2. In this example, the 
mibecraical ly largest component among Zee 240 is Zz, and the largest 
component among ee 256 1s 243: The fault identification region 
R-1 indicates that Eal in NW-1 is deviated less than -20 percent from 
the nominal value. The fault identification region R-13 represents R-] 
in NW-2 because NW-1 has 12 fault identification regions and indicates 
that Eb] is deviated less than -28 percent from the nominal value. Two 
more examples are presented in Fig. 5-4. 

A program for C.E.F.1I., designated by program 2, is written as shown 


in Appendix A and its implementation is annexed in Appendix B. 
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Na oe meee RIMENTS 


A. PASSIVE LINEAR NETWORK EXPERIMENT 
The network shown in Fig. 6-1 (identical to that of Fig. 4-1la) was 


used in the experiment. 


L 


lo /\ L350). Z 
1! , | 9! 


R, = 2 ohms 
L = 2 henry 
C = 3 farads 
Ro = 4 ohms 


Fig. 6-1. 


Passive network 


The following relationships are used for frequency and magnitude 


scaling: 


R* = bDR 
Lees b/ael 
C* = IV/ab C 
f* = af 


where "a" is the frequency scale factor (dimensionless, real and 
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positive constant), "b" is the magnitude scale factor (dimensionless, 
real and positive constant) and a scaled nominal element value is 
identified by an asterisk superscript. 


With a= 10000 and b = 100 , the scaled nominal element values 


were 
1/Eal = R, = 200 ohms 
l/Ea2 = L = 20 milli henry 
Ea3 = C = £3 micro farads 
1/Ea4 = Ro = 400 ohms 


and a set of 3 selected test frequencies became 
Ww) = 20000 radians/sec, f, = 3185 Hz 


Wo = 10200 radians/sec, f, = 1624 Hz 


= 200 radians/sec, f Smo WZ 


The experimental results for fault eermeéions at a set of selected 
frequencies were tabulated in Table 6-1] for the network shown in Fig. 
6-1. Comparing Table 6-1 with Table 4-la and taking into consideration 
the accuracy of the measuring devices, it is concluded that the experi- 
mental results in Table 6-1 and the theoretical results in Table 4-la 
are close. 

Consider the following four cases. 

Case 1. If the resistor, R, » was changed to 130 ohms from the 
nominal value of 200 ohms, the magnitudes of T1l2 were 0.0935 , 0.222 


and 0.755 at a set of 3 selected test frequencies fy ; f, and f. 
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respectively. Since the nominal magnitudes of Tl2 were the same as used 


in Section IV-C, the normalized magnitudes of T12 became 


1.562 


aes 
Xo = e515 
X2 = 1.138 


Case 2. If the inductor, L, was chariged to 30 milli henry from the 
nominal value of 20 milli henry, the magnitudes of Tl2 at a set of 3 
selected frequencies were 0.0462 , 0.131 and 0.664. The normalized 


magnitudes of Tl2 were 


xX) = Oe 7 
Xy = 0.89] 
X2 = 1.0 


Case 3. If the capacitor, C, was changed to 4 micro farads, the 
magnitudes of Tl2 were 0.448 , 0.109 and 0.664 and the normalized 


Magnitudes of Tl2 were 


xX, = 0.748 
Xy = 0.743 
X3_ = 0.98 


Case 4. If it was made that the resistor, Ro» changed to 250 ohms, 


the magnitudes of T12 were 0.0453 , 0.13 and 0.56 and the normalized 


magnitudes of Tl2 were 
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Xy oe Oy 
Xo = 0.889 
X23 = 0.84 


Using the fault identification matrix, Fig. 4-4, Z and the fault 
identification were generated as shwon in Fig. 6-1. 

From Fig. 6-1 and the defined fault identification regions, Table 4-3a, 
the following facts were known for each case. 

Case |. The component Z5 was the largest component in Z and X 
belonged to the fault identification region R-2. Eal was deviated more 
than 20 percent from the nominal value. Ry was out of tolerance. 

Case 2. The network was operating in R-3. Ea2 was deviated to 
somewhere between -30 percent and -50 percent from the nominal value. 

L was out of tolerance. 

Case 3. The network was operating in R-7 and Ea3 was deviated more 
than 25 percent. C was out of tolerance. 

Case 4. Ea4 was deviated more than 50 percent which was indicated 


by R-10. 


B. ACTIVE LINEAR NETWORK EXPERIMENT 

The network shown in Fig. 6-2 was used for active linear network 
experiment. The frequency response of the active network was drawn in 
Fig. 6-3 using the computer program ECAP. A set of selected test fre- 


quencies was chosen as follows. 


Fy = Cee 
f. = 100 Hz 
f, = 1000 Hz 
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100 K ohms 





RY 

Ry = 10 K ohms 

R, =. 1 K ohms 

sm = 5 K ohms 

C, = Syir 

C, = 50 wf 

Q] --- 2N1414 

He = 1300 ohms 
Wa = 36 K aoe 

le = 2.9 10 

He, = 42 

Rea Zed 

GH = Rel" 





(b) small signal ECAP 
equivalent circuit for Q] 


Fig. 6-2. 


Active network 


Since there was no computer program available for fault simulations 
of the active case, the necessary data were obtained from the experi- 
mental results. The fault simulated data were tabulated in Table 6-2. 
The fault simulation curves were drawn in Graph 6-1, Graph 6-2 and Graph 
6-3. The normal regions were defined in Table 6-3 and the fault identi- 


fication regions were defined in Table 6-4. The fault identification 


1 


ee 





matrix was generated by the program 1 (in Appendix A) and is shown in 
fac. 6-5. 

Z and the fault identification were obtained for the following 
three cases. 

Case 1. Let the resistor, Ry » be changed to 14 K ohms from the 
nominal value of 10 K ohms. The measured magnitudes of T12 were 11.71 , 


91.0 and 179.0 at 3 test frequencies f, : fF, and f. respectively. 


Then the normalized magnitudes of T12 became 


Xy = 1.03 
Xy = 06 
(X3 = 1.26 


Case 2. If the resistor, Ro» was changed to 7 K ohms from the 
nominal value of 5 K ohms, the measured magnitudes of T12 were 15.9 , 
126.0 and 190.2 at 3 corresponding test frequencies. The normalized 


magnitudes of T1l2 were 


xy = 1.39 
Xo = ls 
X23 = ioe} 


Case 3. If the capacitor, Ch» was changed to 35 micro farads from 
the nominal value of 50 micro farads, the measured magnitudes of T12 
were 9.5 , 75.0 and 143.0 at 3 corresponding test frequencies. The 


normalized magnitudes of T12 were 
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xX, = 0.834 
Xy = 0.892 
X32 = oOOr, 


Using M, shown in Fig. 6-5, Z and the fault identification generated 
by the computer (program 1) are shwon in Fig. 6-6. 

From the defined fault identification regions, Table 6-4, the follow- 
ing facts were known for each case. 

Case 1. The component Z2 was the algebraically largest component in 
Z. R-3 indicated that R, was deviated more than 30 percent. 

Case 2. The fault identification region R-8 indicated that R. Was 
deviated more than 20 percent. 

Case 3. The fault identification region R-13 indicated that C, Was 


deviated between 40 percent and 70 percent. 


C. PASSIVE AND ACTIVE NETWORK C.E.F.I. EXPERIMENT 

Suppose that the networks shown in Fig. 6-1 and Fig. 6-2 were 
considered as an interconnection of two-port network consisting of two 
cascade connected subnetworks and regarding Fig. 6-1 as NW-1 and Fig. 6-2 
as NW-2, The fault identification matrix for REE experiment were 
Shown in Fig. 6-7. 

Applying the same procedures as in the C.E.F.I. method previously 
discussed, the network malfunctions could be isolated. The C.E.F.I. was 
accomplished for three cases with the same X used in the passive and 
active network single element fault identification experiments. The 


resulted fault identification was shown in Fig. 6-8. 
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Experimental fault simulated data for Active network. 
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DEFINED BOUNDARY POINT VALUES 
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Fig. 6-4. 


Fault identification regions 


printed by computer for Active network. 
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Vili renstoleity OF FAULT PREDICTION 


The feasibility study of fault prediction is quite complicated. 
However a brief discussion of this topic is included in this section. 


Suppose that M of a certain system is 


ae nese 13 an ayy 
Tole 22a 23 27 oor 
avieene. “no. iM” ahr 
Ay Ao A3 Seema = Ss A. Pat inet tea ata iy 
and X' js 
oe I use eX Xx, 1. J 
Z' becomes 
Z' = X'M = [ youn = x, ieee 
= [ Zp 2p mm Be mom, ] 
where 
eee iy 225° aa * j=1,2,--- 47 


x, is positive real number ( i=1,2,---,n) and as; and A. are real 


constants. 
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Since 2 is linear functions of XS; the variation of the value of 
Z. 1S monotonic. 

If the system is faulty and 2. is algebraically the largest component, 
the system is operating in R-j and an element which is out of tolerance 
can be isolated from the predefined fault identification regions. 

Taking the network shown in Fig. 4-la as an example for the fault 
prediction, there are 8 boundary point values between the normal region 
and the faulty region at a test frequency. The largest components in Z 


at these boundary point values are shown in Table 7-1. From this table 


Element xy Xo Xo The_largest component 
(percentage) eraeenndan | in Z at boundary 
Ea) (-20) (0.80 0.80 0.91) eo 
Fal (+20 ) (1.20 1.20 1.06) iee6 
Ea2 (-30) (OneOeeas5 1.02) 1.48 
Fa2 (+25) (1.10 1.05 0.98) 1.28 
Ea3 (-18) ace, lec 1.02) 1.78 
Ea3 ($25) (0.80 0.80 0.98) 1.78 
Ead4 (-20) Viel ele05e 06) 1236 
Ead4 (ciao) (0.90 0.95 0.92) 1.40 

Table 7-1. 


the following observations are made: 

1) If the deviation of Eal element is to be decreased less than 
-20 percent, the largest component in Z must be algebraically greater 
than 1.90; if the deviation of Eal element increased more than +20 
percent, the largest component in Z must be greater than 1.86. The same 


idea can be applied to all elements. 
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Zeke Z for all j is smaller than 1.28, the smallest value in 
Table 7-1, the network is working in normal condition. 

Bait Zs for all j is larger than 1.90, the largest value in Table 
7-1, the network is working in one of the faulty regions. 

4) If the largest component Z. in Z is in between 1.28 and 1.90, 
the network may or may not be in a normal region. For instance, if the 
largest component Z is 1.50 and the network is still in normal condition, 
then it is known that Ea2 and Ea4 are not faulty and Eal or Ea3 may begin 
to approach one of the faulty regions. If the largest component Z. 1s 
1.50 and the network is in faulty condition, then it is known that Ea2 


is faulty. 
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VIII. CONCLUSION 


A fault identification technique has been investigated in this paper. 
Although the examples used are simple, the fault identification proce- 
dures presented here may be extended to more complex systems. 

The theoretical basis of the proposed method is described and the 
technique has been applied to both passive and active networks. 

The main differences in the approach used in this.paper for fault 
identification over the gain signature approach cited in References [3], 
[8] and [10] are summarized as follows. 

1). Fault identification regions are defined uniquely according to 
the fault simulation curves which are constructed from the system sensi- 
tivity instead of assignment of gain signature. 

2). Direct application of the measured values for fault isolation 
can be accomplished instead of gain signature conversion. 

3). Vector-Z indicates an exact faulty element automatically rather 
than a gain signature comparison. 

4). The selection of a set of test frequencies is not critical as 
in the gain signature approach. 

5). By applying the matrix operation to the fault identification 
matrix of the single element fault identification method, one can detect 
not only a fault of the one-element type but also a multi-element fault 
in the given system. 

The proposed method uses only one network function. However, it 


gives no ambiguities as demonstrated in the examples. 
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APPENDIX B 


IMPLEMENTATION OF THE COMPUTER PROGRAMS 


INPUT DATA 
PROGRAM 1 
Card Card 
number column 
1 ] -- 10 
11 -=- 20 
Z ] -- 10 
11 -- 20 
3 ] -- 10 
11 -- 20 
NV+] ] -- 10 
11 -- 20 
NY+2 ] -- 10 
11 -- 20 
2NYV+1 1 -- 10 
11 -- 20 


Information Form 


NV number of selected frequencies integer 


NR number of fault identification 
regions 


lower BP.V. of Xy for fault decima | 
identification region R-1 


upper B.P.V. of x, for fault 
identification region R-1 


lower X9 
upper il Xo il i 
ii il il 
X 
lower - 
u er i xX it TH 
Dp ny 
lower ’ x : ’ 


] 
identification region R-2 
upper i x] 0 i 
identification region R-2 


i il 
lower ay 


upper - 
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2NVt2 


NR NV+] 
NR°NV+2 
NR-NV+3 


NR°NV+4 


NR°NV+2+NM 


] 


LL 


11 
2] 


11 


NV 


20 


10 
10 
20 
30 


10 
20 


lower : XY 
identification region R-3 


upper X 
identification region R-3 


NM number of measurements 


magnitude value of xy 


i 2) 
ul X3 


le 


magnitude value of x 


2 


22 


integer 


decima] 


« 
ul 





PROGRAM 2 


Card 
number 


NR+] 
NR+2 
NR+3 


NR+4 


NR+2+NM 


Card 
column 


2) 


1 


i 


NY 


10 
20 
10 
20 
30 


10 
20 


10 


20 


10 
20 


Infornnation Form 


NI number of M rows integer 
NR number of M columns " 


entry of col-1l, row-1] of M decimal 


row-2 
it row-3 tt 11 
row-ni" i 


entry of col-2, row-1 of M 7 


i row-2 it il 
row-ni : 
NM number of measurements integer 
magnitude value of x; decimal 
il X9 | tt 


ie 


magnitude value of x . 


u 
2 


sy 


* Qpen end value may be used +100. for upper boundary and 0. for 
lower boundary. 
* Appropriate modification may be done. 
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